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. . , Abstract 

, We investigate the possibility of observing very small amplitude low frequency solar oscillations with the 

' proposed laser interferometer space antenna LISA. For frequencies below ^ 2 x 10^^ Hz the dominant 

^SJ I contribution is from the near zone time dependent gravitational quadrupole moments associated with 

^ • the normal modes of oscillation. For frequencies ly above ~ 3 x 10^* Hz the dominant contribution is 

, ^ I from gravitational radiation generated by the quadrupole oscillations which is larger than the Newtonian 

^> • signal by a factor ~ (27rri//c)^, where r is the distance to the Sun, and c is the velocity of light. 

, The low order solar quadrupole pressure and gravity oscillation modes have not yet been detected above 

' the solar background by helioseismic velocity and intensity measurements. Our estimates of the ampli- 

tudes needed to give a detectable signal on a LISA type space laser interferometer imply surface velocity 
' amplitudes on the sun of the order of 1 — 10 mm/sec in the frequency range 1 — 5 10^'^Hz. Such surface 

^1 , velocities are below the current sensitivity limits on helioseismic measurements. If modes exist with 

' frequencies and amplitudes in this range they could be detected with a LISA type laser interferometer. 

■ Keywords: gravitational waves: solar oscillations: space experiments: laser interferometers. 

o : 

. 1. Introduction 

r"| ! The proposed ES A/NASA proposed gravitational wave space interferometer LISA (ESA 1998) consists of 

pi i' three spacecraft at the vertices of an equilateral triangle of sides 5 x 10^ metres; the system is maintained 

^ I in this configuration by arranging that the plane of the detectors has an inclination of 60° to the ecliptic 

5^ ' and counter rotates with the same period as it orbits the Sun. Any two arms constitute a Michelson 

[ type interferometer, a mother spacecraft sends a laser beam to the other two satellites where the signal 

■ is coherently transponded back: the interferometer readout is obtained by interfering the incoming signal 
^ I with the outgoing one and comparing the fractional change in phase shift between the two arms. These 

interferometers will be sensitive to the passage of gravitational waves in the frequency range between 
10~^ and 1 Hz, a range currently inaccessible on the ground due to seismic noise. The primary goal of 
^ . these experiments is to detect waves from individual sources (close binary systems, neutron star or black 

" " ' hole coalescence) and any stochastic background due to the superposition of waves emitted by binary 

systems, - and possibly from the early Universe. 

Current studies of the sensitivity of the LISA experiment (ESA 1998) indicate that the instrumental 
noise, due dominantly to the residual uncompensated accelerations, ~ 3 x 10^^*/\/Hz at 10^"* Hz. In 
addition there is likely to be background "confusion noise" from binary systems which has a comparable 
magnitude at lO""* Hz. (Bender and Hils 1997). When due allowance is made for the angle between the 
arms of the detector and the unknown direction of a source of gravitational radiation, the dimensionless 
strain that could be detected at a signal to noise of 5 is estimated to be ^ 7 x 10~^^ with 1 year's 
observations. 

A LISA type interferometer is sensitive to any variations in the gravitational field in the frequency range 
10^^— IHz. Since the Sun is known to be oscillating in normal modes of small amplitude with frequencies 
in this range, the oscillating external gravitational field would be detectable by LISA if the amplitudes are 
large enough (cf Cutler and Lindblom 1996, Polnarev et al 1996, Roxburgh et al 1999). Only quadrupole 
modes are likely to be large enough to be detected at 1 a.u. and are the modes studied here. In 
addition to the time dependent external Newtonian gravitational field the time varying gravitational 
quadrupole moments also generate gravitational waves, which could give a detectable signal on a LISA 
type interferometer. At a distance of r = 1 a.u. 2-nvrlc > 1 for > f,. 3 x 10~^ Hz, so we anticipate 
that gravitational radiation will be important at frequencies greater than this. 

The first evidence of surface layer solar oscillations dates back to the work of Leighton et al (1962). Low 
order gobal oscillations, which are the oscillations of interest here, were detected by Claverie et al (1979) 
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and Grcc ct al (1980). as resolved peaks in the power spectrum of a time scries of measurements of the 
Doppler shift of a K and Na hne using the integrated hght from the Sun. As a result of ground based 
and space based observational programmes upwards of lO'' oscillation modes have been identified in the 
frequency range 10^^ — 10^^ Hz.; knowledge of these frequencies has been used to infer the acoustic and 
dynamical structure of the Sun (pressure, density and rotation as a function of radius), placing constraints 
on the physics of the solar interior and on models of solar (and thereby stellar) evolution. 

The (p-modes) that have been detected are concentrated in the solar surface layers so the amplitude of 
the oscillating gravitational field is very small. Modes of lower frequency < 5 x 10~^ Hz have amplitudes 
that are relatively much larger in the core and correspondingly a larger external gravitational field for a 
given surface amplitude; since they are more sensitive to the core, knowledge of their frequencies would 
place much tighter constraints on the structure of the very central solar core. Such modes (g- modes and 
low order p-modes) have not yet been detected above the solar background noise which increases at low 
frequencies and great effort is currently being expended in the search for such modes. 

Here we investigate the possibility of detecting such low frequency quadrupole oscillations with a LISA 
type laser interferometer, including both the Newtonian near zone perturbations and the associated grav- 
itational wave emission in the frequency range 3 x 10~^ — 10~^ Hz. We compare the gravitational signals 
detectable by laser interferometry with the velocity signals detectable by whole disc helioseismolgy, and 
demonstrate that low frequency quadrupole oscillations with surface velocity amplitudes below current 
helioseimic limits could nevertheless be large enough to be detectable by LISA. If such modes are first 
detected by helioseismic techniques then the measured frequencies (and predicted power) would provide 
a valuable calibration tool for LISA; if they have not been detected by helioseismic means then one can 
look upon LISA as a potential telescope for studying the deep solar interior. 

The plan of this paper is as follows: first we discuss the properties of the solar oscillations and the 
quantify the relationship between the surface radial amplitude of an oscillation, its quadrupole moment 
and its horizontal amplitude. We then express the external gravitational field in terms of the quadrupole 
moment tensor expressed as a sum over a set of basis tensors (corresponding to surface harmonics) which 
enables us to relate the generation of gravitational waves to the surface amplitude of the oscillation. We 
then determine the response of a LISA type detector to both the time dependent Newtonian field and to 
the associated gravitational waves for a given oscillating quadrupole moment; above ^ 3 x 10~^ Hz the 
signal is dominated by gravitational waves. Next we determine the magnitude of the velocity signal from 
each mode in terms of surface amplitude. We then consider the background noise for both velocity and 
gravitational detectors and hence the signal to noise in both gravity and velocity for a given assumed 
frequency resolution. The ratio of these S/N is then independent of the assumed surface amplitude of 
the oscillation and if greater than 1 the modes are easier to detect by a gravitational laser interferometer 
than by helioseismic experiments. The outcome of these calculations is summarised in Figure 6, For 
frequencies < 2 x 10~^ Hz the m = 2 (and m = 1) quadrupole modes are more readily detected in by a 
LISA type interferometer than by helioseismic experiments. 

2. Solar Oscillations and space-time around the Sun 

The solar oscillations are normally expressed in terms of a surface harmonic and Fourier time decompo- 
sition with any variable, eg the displacement 6r expressed in the form 



where uj = ujnem are the cyclic cigenfrequencies of modes corresponding to a particular surface harmonic 
S i.m(() 1 (j)) ^ n labels the order of the mode (|r7,| is essentially the overtone number, the number of nodes 
in the radial direction), Cnfm(^) the corresponding dimensionless eigenfunctions and (r, 6*, </>) spherical 
polar coordinates with origin at the centre of the Sun. The modes are classified as p (pressure) modes 
with frequencies increasing with increasing n, and g (gravity) modes with frequencies decreasing with 
increasing n, (for clarity we use negative n for the g-modes in equation (1)). If the basic miperturbed 
state is spherically symmetric the frequencies are independent of azimuthal order m. Rotation lifts this 
degeneracy giving frequencies uonim ~ w„£o + rnO, where fl ~ ~ 3 x 10~^ rad/sec is a weighted mean 
of the solar angular velocity. With a frequency resolution of Aiy ~ 3 x 10~^ Hz, as envisaged for LISA 
(with 1 year's observation time), these individual m- value modes should be resolved provided there is 
sufficient power and the line widths are sufficiently narrow. For frequencies u = co/2it 10~^ Hz, which 
is the region of interest in the present analysis, Qq/oj 10~^ and the eigenfunctions of the modes may 
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be taken to be independent of azimuthal order m. The solar rotation axis is inclined at an angle ^ 7° to 
the ecliptic plane, which introduces an additional (small) modulation which is neglected in the present 
analysis. 

Since oscillations with frequencies ~ 10~^ Hz have not yet been detected in the solar oscillations, the line 
widths are Tinknown; the measured line widths at higher frequencies decrease with decreasing frequency 
and crude extrapolation from the measured range suggests Av < 3 x 10^* Hz at frequencies ^ 10^^ Hz 
and consequently mode lifetimes > 1 year. We shall assume here that the lines are narrower than the 
frequency resolution and the modes may effectively be considered as monochromatic. Only the quadrupole 
modes i — 2 arc considered; the external gravitational potential of a multipole of order i decreases like 
l/r^~^^ and modes with £ = 2 will dominate at 1 a.u. Further since the oscillation velocities are very 
small compared with the velocity of light only quadrupole radiation will be significant. 



The external Newtonian gravitational potential of the oscillating Sun can then be expressed in the equiv- 
alent forms 

U{v, t) = -G j ^ = f/o(r) - ^ V-P (i) = C/o(r) - GM^RI ^ '^S^mie, </>) (2) 

^ ' m=-2 n=-oo 

where C/o(r) is the time independent potential, surface harmonics of degree £ = 2, J^m oc e**^* the 
dimensionless time dependent quadrupole moments corresponding to eigenmodes with cyclical frequencies 

V"^= [ {3x"x^ -6'''^x''x^,)p{t)dV (3) 
Je 

is the quadrupole moment tensor, a Cartesian coordinate system with = the x^ axis and (p = 
the x^ axis, 6"^ is the Kronecker delta and Va/j = d'^/dx°'dx^ . 

We take the as real surface harmonics normalised to unity over the sphere 



5 3 cos^ 9-1 ^ I 15 cos / 15 . ^ n cos 2(j) 



introduce the set of trace-free basis tensors (corresponding to surface harmonics) as 
10 

Jo = I 1 I Ji = I I J_i = I 1 I J2 = I -1 1 2:_2 = I 1 1 (5) 
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and expand the quadrupole moment tensor as 
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V"^ = MqRI J2 H CmJnml^^, where Co = -^ — , Cm = \j^ m^O (6) 

n= — 00 m=— 2 

The relationships between the S2m,Im, and Jnm are derived in Appendix 1. 

The properties of the eigenmodes were computed using the (iwr-smoothed) standard Aarhus Solar Model 
SI of Christensen-Dalsgaard et al (1996), with the dimensionless radial eigenfunction Cn(-R) normalised 
to 1 at the solar surface r = Rq. The amplitude of oscillating quadrupole moments J„ are given in Figure 
1. (We here suppress the subscript £ — 2 and m, since the eigensolutions are independent of m.) Table 
1 gives a summary of the results obtained for the amplitudes of relevant quantities derived in this and 
subsequent sections: column 1 gives the radial order n of the mode, column 2 the frequency v, column 
3 the horizontal displacement eigenfunction at the solar surface (h, column 4 the quadrupole moment 
J. The remaining columns give the computed signals for these modes as evaluated in the following 
sections: column 5,6,7 gives the amplitude of the gravitational signal function Sm for m = 0, ±1, ±2, and 
in columns 8, 9 the amplitude of the velocity signal Vm for m = 0, ±2. Note these are calculated with 
the normalisation Cn(-R©) = 1; Ch, J, Sm, Vm scale linearly with the value of ({R). We observe here that 
since the m = ±1 modes are antisymmetric about the equator plane 9 = Tr/2, these modes will have no 
net velocity so = 0,m = ±1. These modes do however have a gravitational signal. 
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The sun is rotating with a period ^ 27 days so a frame of reference fixed relative to the osciUating sun 
rotates relative to an inertial frame, producing a modulation of the time dependent gravitational field. 
The solar rotation axis is inclined to the ecliptic plane, and hence to the orbit plane of the gravitational 
detector, producing a further modulation of the signals from the oscillations and a very low frequency 
signal (~ 4 X 10~^ Hz) from the static quadrupole moment induced by the rotation. As the inclination 
is small (~ 7°) we neglect these effects in the present analysis taking the rotation axis of the Sun 
perpendicular to the orbit plane of the detector. 

Since the monopolc gravitational field of the Sun is weak {GMQ/rc^ ~ 10~^ at 1 a.u.), and that of the 
time dependent quadrupole moments and any associated gravitational radiation even weaker, the metric 
of space-time in the neighbourhood of the Earth is adequately described by the weak field limit as 

ds^ = Qikdx' dx^ = {rjik + hik)dx' dx^ = + 5") '^^'^^^ ~ ~ 5") ^^"'^^^ + h^fdx'^dx'^ (7) 

where ry^fe = diag{l, —1, —1, —1) is the Minkowski metric, \hik\ << 1, and indeces are raised and lowered 
using the Minkowski tensor. Roman indeces i,k = 0, 1, 2, 3 whereas Greek indeces a,(3 = 1, 2, 3. The hik 
have a contribution from the time dependent quadrupole moments 



[l)=-^^6.kV.Jl) ± tc^Jr^^T-J (8) 

^ ^ ^ ^ n= — oo m=— 2 

and a space-like contribution from the gravitational quadrupole radiation, which with time dependence 



25U _ GMqR% fl\ GM^Rl 

"-ifc ~ 3^2 "^'^ ^ 



~3^d^ V "^J ~ ~~3^ " 3^ ^ ^ C^Jn-mda^he^^m (9) 



n=— oo m=— 2 



where Vap is the retarded, transverse, trace-free part of Vaii (ie in the plane perpendicular to the 
radial direction of propagation of the waves) and are the retarded, transverse, trace-free tensors 

corresponding to the basis tensors . 

In this weak field limit the hik can be considered as a superposition of independent components and 
the effect of each component considered independently. In the following section we drop the indeces n 
which label the radial order of modes and the index m which signifies the basis tensor (or equivalently 
the spherical harmonic). 



3. Phase shift in a perturbed space-time 

Consider an electromagnetic wave of frequency uoe and wave 4- vector ki propagating along the arm of a 
detector from an emmitter at x\ to a receiver at a;^. The phase (f) of the wave at a;^ is given by 



f 

J A 



B 

kidx'' = {kiBXg - kiAx\) - I x^dki = kiA{xB-x]i)+ I {xg - x^)dki (10) 



In Minkowski space the wave vector A:* is constant = k\ = ci;e(l,n")/c and the phase difference cf){B) — 
(j){A) is simply kiA{x^B ~ ^k)- The second term on the right hand side of equation (1) is identically zero 
in the Minkowski limit. 

In the time-dependent space time gik the first order perturbation in phase, due both to the departure 
from Minkowski space time and the displacement of the receiver B (taking A as fixed) is given by 

5ip = kiAS{x% - x\) + [ {x% - x')dki (11) 

J A 

We take A as the origin of coordinates so a;^ = 0, and take the unperturbed ray to be given by a;* = A/c', 
where A is an affine parameter varying from at A, to A at B. Since k^ = dx'^/dX, the null geodesic 
equation in this weak field approximation reduces to 

= — r^^fc^fc" Ri ■^hmn,ik"^k" (12) 
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and the last term in equation (11) can be expressed as 

f-A 1 /-A 1 /-A 



y (A - X)k'dki = -J (A - X)hmn,ik"'k''k'dX = -J (A - X)-^k"'k'^d\ (13) 

Substituting this result into equation (11), defining V = (x^ — a:^), (5/* = (5(a;^ — x\), and integrating 
by parts gives 



5ip = kiASP + 



^{A-X)hmnk"'k'' 



A ^ 



2-^0 



hmnk"'k''dX (14) 



is given by the geodesic deviation equation 

:(<5r)=i?k„ — — ^^ which reduces to _(<5r) = i?Uo«' (15) 



£>t2 ^ ' ^""^ dr dT ' df 

in the slow motion approximation appropriate to the current analysis. Since in the weak field approxi- 
mation 

Kko = Kk,o - Ko,k + ^iiKk - « Iv'' {hki,oo - hoi,ok + Ko,ik + hok,io) (16) 

it follows that = 0, d'^5l°/dt^ = 0. RecalUng that hoa = 0, I' =1^+ Sl^ with I'q = (1, n) we obtain 



1 r^2 



dt^ 



IS (17) 



where again Greek indeces denote spatial components and V^^ = d^/dx'^dx^. Since the time dependence 
of hik oc e"^* this equation can be integrated to give 

SI" = - ^"^i^f^oo) ll (18) 
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Now since fc" = uieu" /c^ A = cl/uje, where We is the frequency of the electromagnetic wave, substitution 
into equation(14) gives 

= ^ (/ia/3n«n'5 - ^n^nNl0ho)j - ^/l„^„n'"n" + /i„„n™n"rft (19) 

and / = ria^o • 

We now write = 5(/5Ar + 5(pGW where 5^pn is the contributions from the Newtonian potential 5U and 
5ipGW that from gravitational waves and expand the hmn in the last term by Taylor series to give 

+luJen'^Vahoo [ ' e'^'^ctdt + -uj^n'^n'^Vlphoo I ' e^^cH^dt + . . . (20) 

S^GW = Icoe r h^^n'^n^dt « \ ^n'^n^h^J (21) 
^ Jo ^ c 

We have here used the fact that in the weak field limit h^^ — hoo^ap (see equation 8) to eliminate ft^^. 
Retaining just the leading terms in powers of col/c« 1 and l/r « 1, the phase shift reduces to 



where V^^^^ = d'^ / dx°'dx^ dx^ dx'^ and we have used the results (8) and (9) for ft,^ and hf^. Since the 
unperturbed phase shift ip = ix>el/c the fractional change in phase shift can be expressed in the form 

where Nai3 = rian^, and 

Ta/3 = ^ r^N'^'Vip^^ = 2nan0 - lOAtKn/j + narip) + 5(7ai^ - (24) 

where is the unit vector in the radial direction and fi = n^na- The details of the derivation of T^p 
are given in Appendix 2. 



4. Response of the interferometer 

A laser interferometer detector such as LISA consists of 3 arms AB, AC, CB, in circular orbit aroimd 
the Sun; the interferometer's response is given by the difference in the fractional change in phase shifts 
of the round trip signals along any two of the arms eg. ABA and ACA, as 



AB 



G 



AC 



2wV5 



(25) 



where A indicates the difference between the arms AB and AC. 

The response s is a function of the angles x, l, tp, v, the quadrupole moments Jm and the frequencu u. x 
is the angle between the radius vector to the detector and the instantaneous orientation of the reference 
axes corotating with the sun, t the inclinations of the plane of the detector to the orbit plane, tp the 
angle between the arm AB and the direction of the orbit and v the angle between the two arms of the 
detector. In the LISA concept i = v = n/S and the detector rotates with the same period as it orbits 
the Sun, ip decreasing as cj) increases, where cj) is the position angle of the detector on its orbit relative to 
a fixed inertial frame. 

We now express the quadrupole moment tensor in terms of the trace-free basis tensors (see equation 
6) and determine the response Sm for each mode of azimuthal order m and frequency v = u)/2'jt 



1 fRc 

2 \ r 



GW 



where 



al3 




f— ) 
\2iTr) 



3 X IQ-^Hz 



(26) 



and are the transverse-trace- free components of I';^ ■ We note that provided and are of 

comparable order the gravitational wave contribution dominates for v > ^ 'i x lO^^Hz. 

To determine the functions consider the detector to be at the point P on its orbit with orbit 

angle (j) relative to a fixed inertial frame, and let x the angle between the radius vector to P and the 
axis of the corotating system x". We then define local transverse Cartesian coordinates ^'^ at P such 
that is in the outward radial direction, in the direction of the orbit and perpendicular to the 
orbit plane. The basis tensors in the ^'^ coordinates are given by 



where = (cosx,sinx,0), = (- sinx, cosx, 0), e^ = (0,0,1) 



(27) 



are the unit vectors along the axes in the coordinate system. The transverse trace-free tensors T^^ 
are then given by (sec Appendix 3) 



t23 tJ^ T. 



32 



t23 



^22 _ 



_t33 _ 1 /t22 _ t33\ 
m 2 



= 



(28) 
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The detector lies in a plane inclined at an angle l to the orbit plane, the arm AB making an angle tj) to 
the orbital direction and AC an angle ^ + u. In the ^'^ coordinates the unit vectors n^, are then: 



(— sin ^ cost, cosV'jSin'^sin t), Uq = (— sin^' cost, cosip', sin^' sint), ip' = ip + v 



(29) 



The unit radial vector fi"^ = (1, 0, 0) and after some laborious algebra we determine the tensors AN^^, ATj^ 
in the ^'^ coordinates, and hence the functions f^^fm^ defined in equation (26). The solutions for 
arbitrary L,v,ipjX ^re given in Appendix 3. 

In the LISA concept l = v = w/S and the detector rotates with the same period as the orbit in the sense 

that ip decreases as (/) increases, so that ijj ^ tpQ — (f) where iJjq is a constant, the orientation of the detector 
relative to the orbit plane at the arbitrary zero of the orbit angle 4>. We set tpo = and obtain 



Jo 



N 



3\/3 



sin(2^ + 7r/3), 



GW 



2lV3 
16 



sin(2V' + 7r/3) 



3 

/j^ = -3 [ 2 cos X sin(2V' + 7r/3) - sin x cos(2V' + 7r/3) ] , /f ^ = - sin x cos(2V' + 7r/3) 



rN _ 



25 

— cos2x sin(2'^ + 7r/3) — 8 sin2x cos(2'^ + n/3) 



fGW 
/2 



7V3 
1^ 



cos2x sin(2V' + 7r/3) 



(30) 



The values for m, = —1, —2 are obtained from those for m = 1, 2 by replacing x by x ^ jlm. 
Relative to a fixed inertial frane at time t the detector is at orbit angle </> = and the corotating 
axis is at an angle fi©*, where il^ is the angular velocity of the detector around the sun and the 
angular velocity of the Sun relative to this inertial frame. The angle between the radius vector from the 
sun to the detector and the rotating axis is therefore x = where = 2it/Ps where Pg ~ 26.75 
days is the synodic period of solar rotation, that is the period relative to a reference frame orbiting the 
sun at 1 a.u. Substituting into the results (30) gives /„((/') and fm^i'f') a-s a function of the position of 
the detector on its orbit determined by the orbit phase angle </>. The f^W are shown in Figures 2a-2c 
and the f^^{(f>) in figures 3a-3c. Since the m = mode is axially symmetric the /o(^) are independent 
of the rotation of the Sun, the other modes display the modulation of the signal due to solar rotation. 
Combining the contributions from and gives 



So{l',(j)) = -T^CqJo — 

16 r 



7 z/ 
3 



sin(2V' + 7r/3) 



1 

1 + - — ) sinx cos(2V' + n/3) - 2 cos x sin(2V' + tt/3) 



R5 



sin2x cos(2V' + 7r/3) - ^ ( 1 + ^ ^ ) cos2x sin(2V' + 7r/3) 
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where tp = -(t>, x = (t>Pd/Pe 



13.65 ( 



Co 





m ■■ 



±1,±2 



(31) 



(32) 



Again values for m = —1,-2 are obtained from those for m = 1,2 by replacing x by x ~ 7r/2m. The 
results (31) give explicitly the contributions from the Newtonian potential and gravitational waves and, 
as anticipated above, the gravitational wave contribution dominates when v > i/r ^ 3 x 10~*IIz. 

The amplitude of the signal of given (m, v) varies over the course of a solar month and during the year. 
For an observation time T ~ 1 year the amplitude is given by the root mean square averaged over T. The 
resulting amplitudes are given in Table 1 (columns 5,6,7) for values of the quadrupole moments given in 
column 3 corresponding to unit surface radial eigenfunction. These signals scale linearly with the value 
of the quadrupole moments. 



5. Surface velocity amplitudes on the Sun 

The solar oscillations of low degree are obtained by taking the power spectrum of a time series of mea- 
surements of the Doppler shift (or velocity) of a given spectral line in the integrated light from the sun 
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- the frequencies are the peaks in this power spectrnrn. The velocity amphtude of a given quadrupole 
oscillation is therefore given by integrating the component of the surface velocity in the direction of the 
observer over the visible solar disc, 

_ Jk°'Vafin)ndS 

where in this section we define fc" as the unit vector in the direction of the observer, the cosine 

of the angle between the unit radius vector fa and the direction to the observer, and /(yu) the appropriate 
solar limb darkening function which incorporates the angular dependence of the intensity of radiation at 
the solar surface (see for example Allen 1973). We here take /(/it) = a + bn with a = 0.55, b = 0.45, which 
is a reasonable approximation for the whole disc velocity measurements by the GOLF (Turk-Chieze 1998) 
and BiSON (Chaplin et al 1998) experiments 

In a system of spherical coordinates (r, 6, (f>') corotating with the sun, the velocity on the solar surface 
due to an oscillation mode is 

/ d 1 d \ 

Va = iVr,Ve,V4,) = (J Rq Ur, (hg^, Oi^^q^ j S2m{0,(p') (34) 

where (r, (h are the dimensionless radial and horizontal displacement eigenfunction at the solar surface 
R = Rq . Let the corotating axes be at an angle 0o to the direction to the observer and take (j) — (p' — (po 
then 

k"va = Vr sin cos (j) + vg cos 6 cos (f) — sin <p, /X = sin ^? cos (/> (35) 
and evaluating the integral (33) gives the velocity amplitudes 



Vo = -i 



fV / (16a+15b)C. + 6(8a + 5b)a \ /;- cos2.^o _ , o^ /o7^ 

W^l 40(3a + 26) J' ^±^ = ^^^° sin 2,^0 " ^^^^ 

For m = ±l,52m oc cos^ is antisymmetric about 6 = 7r/2 and the integrals (32) for V±i are identically 
zero. The m = ±1 modes are not detectable in integrated velocity in the equatorial plane. 
The surface horizontal displacement (.h is known in terms of the surface value of the radial displacement 
eigenfunction (^r from integration of the eigenvalue equations for the oscillation (column 3 of Table 1). 
Since the normalised value of the quadrupole moment (column 4) is << 1 for all modes, it follows from the 
equations governing the oscillations that (h w C.Rg/uj'^ Rq (see Unno et al 1989). The velocity amplitudes 
with (r = 1 are listed in column 3 of Table 1; they scale linearly with the value of C,r. 



6. Background Noise the Gravitational and Helioseismology Experiments 

The cxpcictcid sensitivity of the LISA experiment for a single 2-arm detector has recently been rc!-c;valuated; 
at low frequencies the dominant contribution is from the acceleration noise oc 1/w^, at intermediate 
frequencies from shot noise, and at high frequencies the sensitivity declines when the path length along 
the detector becomes comparable to and greater than the wavelength of the gravitational wave (ESA 
1998). In the region of interest this background noise level is given by 

/ 1X^/2, 

Si(i/) = ( 1.6 X lO"'^^ + 9.2 X 10"^^^ j /VW^ (38 

This background noise is normally presented in terms of limits on a signal that could be detected with 
one years data at a sampling rate of 1 sec, with a signal to noise of 5, for an unknown location of the 
source, and including the inclination between the two arms of the detector. To reproduce the limits given 
in ESA (1998) we divide the result (38) by VT (in sees), mutiply by 5 for the S/N, by \/E for the unknown 
direction and divide by sin7r/3 for the angle between the two arms. This is shown in Figure 4. 

In addition to this instrumental noise there is also a contribution from galactic and extragalactic binary 
systems, the magnitude of this "confusion noise", Bi,{v), is uncertain but has been estimated by Bender 
and Hils (1997); their estimate is given in column 11 of Table 1 and shown in Figure 4. 
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The background noise in velocity experinw^nts is doininantly from velocity fields on the solar surface (active 
regions, granulation, meso-granulation, supergranulation) and the cumulative effect of these motions has 
been estimated by Harvey (1995). However in contrast to the gravitational case this background noise 
has been measured in a number of helioseismology experiments and in Figure 5 we show the background 
velocity noise determined by the GOLF experiment on SOHO (Gabriel 1997,Turck-Chieze 1998). This 
is in reasonable agreement with the value obtained by the ground based whole disc networks BiSON 
(Elseworth et al 1994) - who obtained 900 cm/sec/\/Fz at v = 5 x 10"'' Hz and IRIS (Fossat et al 1996) 
who obtained 2 x 10^ cm/sec/ VWz a,t u = 10~^ Hz. These values are a factor ~ 5 below the model 
predictions of Harvey. In the region of interest this backkground noise can be approximated by 



7. Comparison of detectibility in gravitational and velocity experiments 

We now make the following comparison. We assume a frequency resolution in the power series analysis 
of 3 X 10~^ Hz (corresponding to 1 year's observation time) and the value of a velocity signal from an 
oscillation mode that has S/N=3 using the above background noise estimate (39). This gives a value 
of the dimensionless radial eigenfunction through equations (37) for the velocity amplitude for each 
frequency and m = 0,2. This then determines the amplitude of the quadrupole moment J and the 
gravitational signal strength S obtained by simply scaling Sm in Table 1 by the value of Ci?- We then 
compare this gravitational signal strength with the background noise for the LISA experiment computed 
from (38) under the same assumptions: Azv = 3 x 10~^ Hz, S/N=3. The results are shown in Figure 
6. We see that for frequencies v < 2 x 10~^ Hz the signal from the m=2 modes are stronger in the 
gravitational experiment than in the velocity experiment On this diagram we also show the S/N for the 
gravitational experiment when binary confusion noise is included. Note that since the averaged velocity 
for the m = ±1 modes is zero there is no limit on the gravitational signal set by the helioseismology 
experiments. 

This comparison between the 2 experiments is of course independent of the assumed S/N in the velocity 
experiment, what is being compared is the detectibility of an oscillation mode by the two techniques and 
is essentially the ratio of the S/N for the gravitational experiment to S/N for the velocity experiment. 

We can reverse the comparison and ask what would be the velocity amplitude of modes that are at 
the margin of detectibility with the gravitational detector, and how does this compare with the solar 
background noise? This is done in Figure 7, here we see that the velocity amplitude decreases at low 
frequencies and for the m = 2 modes is well below the noise level for a frequency resolution of 30 nHz. 
The local minimum at ~ 4 x 10~'^Hz is due to the inclusion of gravitational wave emission. 

In the above we have assumed that the signal is monochromatic - or rather that the line is narrower than 

the frequency resolution. This is the prediction from extrapolating the p-modes observed by helioseis- 
mology at higher frequencies. Since the rotational splitting of the modes is known, signal enhancment 
techniques can be used, superposing power in frequency bins separated by mfi©, as is being done in the 
search for g- modes in helioseismology by the Phoebus group (Frolich et al 1998). If the modes have a line 
width in excess of the frequency resolution the detectibility in both velocity and gravity is correspondingly 
reduced but the ratio of S/N for the two experiments remains the same. 

8. Conclusions 

At low frequencies, v < 2.2 x lO^^'Hz, the quadrupole oscillations of the Sun are more readily detected 
through their gravitational signals on a LISA type space interferometer than by helioseismology. The 
emission of gravitational waves by the oscillations contributes to the gravitational signal for frequencies 
> 3 X 10-4 Hz. 

These low frequency modes are more sensitive to the structure of the solar core than the higher frequency 
p-modes that have so far been measured by helioseismology experiments. Were these low frequency modes 
detected they would advance our understanding of the structure of this cxnitral core. The helioseismology 
experiments measure the surface velocity which is dependent on the detailed structure of the outer layers 
of the sun, and this has to be subtraced off to give diagnostic information on the solar interior. The 
interferometer experiments however measure the actual quadrupole moments which are determined by 




(39) 
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the oscillations in the deep high density solar interior and therefore encode information about the solar 
interior which is independent of the structure and physics of the outer layers. 

The predicted amplitudes of these low n modes is also uncertain. The p-modes are thought to be 
stochastically excited by the convective motions, and if this is the excitation mechanism of the g-modes 
the predicted amplitudes of the low frequency modes are very small. But as was first pointed out by 
Dilke and Gough (1972) and Christensen-Dalsgaard et al (1974), the steep gradient of ^He in the inner 
half of the sun could also provide an excitation mechanism for these modes which may be damped by 
parametric resonance with other modes or by mild turbulent diffusion, but still be of sufficient amplitude 
to give rise to a detectable gravitational signal. 

We note also the the frequency of 10"''' Hz corresponds to a period of the order of 160 minutes. There 
have been repeated suggestions that such a signal has been seen in helioseismology experiments, many 
but not all of the claims being withdrawn following more detailed analysis (but see Kotov et al 1997). In 
a recent analysis of the GOLF data Palle et al (1998) place an upper limit on the velocity amplitude of 
such a mode ~ 1 cm/sec, an estimate compatible with searches by the Phoebus group (Frolich et al 1998). 
We note that for such a 160 min mode an S/N=l in the helioseismic velocity experiments corresponds to 
a S/N=20 in a LISA type experiment. 

The search for g-modes in the helioseismology experiments is ongoing and hopefully will detect some 
modes prior to the launch of LISA. Were such modes to have been identified then they would provide 
a valuable known signal on a spac;e interferometer which could then be used to calibrate and test the 
experiment, giving more credence to the interpretation of signals from more distant astrophysical sources. 
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Appendix 1. The quadrupole basis tensors X"^ and surface harmonics S2m 

The quadrupole gravitational potential can be expressed in the two equivalent forms 

GM^RlY, E ^^2„^(0,</>)=Gi2?"^V„^(i) (^1) 

m=— 2 n=— oo ^ ^ 



where (r, 0, (p) are spherical polar coordinates, = {x, y, z) Cartesian coordinates with 9 = the z axis 
and (j) = Q the x axis, and Va/3 = d"^ /dx°'dx^ . Jnm are the (dimensionless) quadrupole moments, X>"^ the 
quadrupole moment tensor and S2m are the real surface harmonics normalised to unity over the sphere 



„ /5 3cos^^-l / 15 . „^cos(2!) „ / 15 . 2/,cos2c!) , 

V 4-77 2 ' V 167r sm0' ' V 167r sm20 ^ ^ 

We define a set of 5 independent trace free tensors as 

/I \ /O l\ /O 0\ /I o\ /o 1 o\ 

Jo = 1 Xi = J_i = 1 J2 = -1 J_2 = 1 (A3) 
\0 -2/ \1 0/ \0 1 0/ \0 0/ \0 0/ 

Now since 

l\ 9^ /1\ f SXaXff ~ r^da/n 



where 5c«/3 is the Kronecker delta, and 

= X = rsinO cos 0, x^ = y = rsinO sin 0, = 2 = r cos 6, {A5) 

the satisfy the relations 

^ /1\ 3,9 , „n, 6 1 — 3cos^^ 6 i4jt „ , 
If W.A-]=-^{x^^y^- 2z') = -3 = --,J^ 52,0 (A6) 



ir'^Va/j ( - ) = ^22;^= -1 sin^cos^cos.^= 4a/^52,i {A7) 



r r 



2:"?Va/3 ( - ) = ;::5 2y^= ^ sin^cos^sin</.= -3 y^52,-i (A8) 



^^'^ v., { i) = 1 (x^ - y^) = A sin^ ^ cos 2<^ = ^ ^/ ^ 52,2 (A9) 



^-^ (i) = I 2xy = A W ^sin2<^ = ^ V ^ (^10) 
If we now expand the quadrupole moment tensor as 

oc 2 

V^I' = MqRI J2 C'mJnml^^ (All) 

n= — OQ m— — 2 

then on substituting into eqn (Al) the coefficients Cm are found to be 

^o = -\/^' C'm = \/^ m = ±l,±2 (A12) 

V 477 V 477 
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Appendix 2. Derivation of the detector tensor T^^ 

We use a cartesian coordinate system a;" with x'^Xo. = r^- Then 

v.ll 



V,,, ll] = +3^ + 3^ + 3^ - 15^« 



1\ 3 



Defining the unit vector in the radial direction 

Xi X2 X3 

rie = ( — , — , — ) 
r r r 

then 



Defining = n'^n^ and contracting with n'n''' gives 



Since this tensor is to be contracted with the symmetric trace free tensor X"^" = I?"*^ and SfSaT)' 
we only need the trace fee component of this so the terms with Sa/3 can be removed leaving 

Ta/3 = Y'^'^'^'^'^afj-yS f^) = 2n/3na - lO/OT/sna - lO/inan/3 + 5(7/U^ - l)n/3na 
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Appendix 3. The functions /f (^), ff^{(l>). 

Here we compute the source tensor in local transverse Cartesian coordinates ^" at a point P on the 
orbit of the detector where are in the orbit plane, in the the outward radial direction, is in 

the direction of motion and perpendicular to the orbit plane, x is the angle between and the 
direction of the coordinate system in which the multipole moments are determined. The unit vectors 
along the ^" axes in the x"' coordinate system are 

4 = (cosx, sinx, 0), = (- sinx, cosx, 0), e| = (0''^'l) 

Hence the basis quadrupole tensors in the ^ coordinate system are given by = e^^e^I"'' and are 
















cosx 


-sinx 










0/ 








cos 2x — sin 2x 
I2 = \ sin 2x — cos 2x 

oy \ 

The transverse-trace-free radiation tensors in the ^ coordinates are then 





smx U U 11 i_i =cosx 




^ /O \ ^ /O 

l2 = -^cos2x 1 L2 = --sin2x 1 
^ \0 -1/ \0 -1 

Now let the plane of the detector be at an inclination angle l to the plane of the orbit (t = tt/3 in the 
LISA experiment), and let the detector arm AB be at an angle to the direction of the orbit (^^). In 
the the ^ coordinate system the the unit vector along the arm of the arm AB of the detector is 

n" = {—sintp cos t, cos tpjSintp sin t) 

so the projection tensor iV^^ = n°^n^ is 

sin^ ij} cos^ L — sin ip cos ip cos i — sin^ ^p sin l cos l 
Nab = I — sin cos tp cos l cos^ ip cos tp sin ip sin i 

■ sin^ Ip sin t cos t cos ip sin sin l sin^ V sin^ i- 

Lct the angle between the arms AC and AB of the detector be v so that AC makes an angle tp + v with 
the direction, the projection tensor N^c is given by the above result with tp replaced by + w so the 
difference 

(— sm{2tp + v) cos^ L cos(2i/; -|- v) cos t sm{2tp + v) sin l cos l ^ 
cos(2i/' + u) cos t sm{2tp + v) — cos(2-!/' + w) sin /, 

sin(2i/' + v) sin l cos i — cos(2'0 + v) sin l — sin(2'0 + v) siia^ t 

The signal functions for gravitational waves are given by 

f^^ = ANc^il^^ = sin V (sm{2ip + v) 1^ - 2 cos(2V' + v) sin l - sm{2ip + v) sin^ 1 1^^ 



which g 


;ives 






fGW 
JO — 


1.5 sinu (1 + sin^ l) sin(2-!/) 


+ v) 




fGW 
Jl — 


2 sin V sin l sin x cos(2V' + 


v) 




fGW 
J-1 — 


—2 sinw sini cosx cos(2V-' + 


v) 




fGW 
J2 — 


—0.5 sin V [1 + siv? l) cos 2x 


sin(2V' - 




fGW 
J -2 — 


—0.5 sin v (1 + sin^ l) sin 2x 


sin(2i/) -\ 


-v)) 
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The Newtonian signal functions = ATaffl^^ where 

AT„/3 = 2AN„p + 10AM^i3 + 5Ai?a/3 

with 

Ma0 = fi{n°'n'^ + fi'^n"), Ra0 = {Iji^ - l)n"n^, n = ffn-y 

In the ^ coordinates n" = (1, 0, 0), and for the detector arm AB ji = n"na = — sinipcosL so 

2 sin^ ip cos^ i — sin i/) cos ip cos l — sin^ ^ sin l cos l ^ 

- sin tl) cos tl) cos i 

- sin^ ip sin i cos l 



— 2 sin(2?/' + w) cos^ i cos(2'^ + w) cos (. sin(2'^ + i;) sin t cos (- ^ 

cos(2i/' + cos (. 

sin(2'^ + v) sin t cos t 

The tensor i?^^ = {T/j."^ — l)n"n'^ gives 

/ 1 0" 

Aii = -7sint;sin(2V' + v)cos^i 

\0 

So finally the tensor AT^^ is 

(— 17sin(2V' + u)cos^(- — 8 cos(2Vj + -u) cos t — 8 sin(2V' + sin tcos 
— 8 cos(2?/' + cos (. 2sin{2^ + v) — 2 cos(2V' + f) sin t 

— 8sin(27/; + sintcosi — 2 cos(2^ + sin t — 2 sin(2i/> + sin^ t 

The signal functions = I^^ATa/s are then 

= sinv (6 — 21 cos^ l) sin(2V' + v) 
fi = — 4 sin t; sin l [4 cos l cos % sm{2il> + v) — sin x cos(2?/; + 1;)] 
f^i = —4 sinu sini [4 cost sinx sin(2?/' + v) + cosx cos(2-0 + v)] 
f2 = — sinw [(17cos^ L + 2) cos2x sin(2V' + v) — 16 cost sin 2% cos(2V' + f)] 
f!^2 = — sinv [(17cos^ t + 2) sin 2% sin(2'^ + v) + 16 cost cos 2% cos(2^ + u)] 
setting L = V = tt/S gives 

/o^ = ^ sin(2^ + 7r/3), = ^ sin(2^ + ^3) 

/f = -3 [2 cosx sin(2V' + 7r/3) - sinx cos(2V' + 7r/3) ] , /f ^ = 1.5 sinx cos(2V' + 7r/3) 
/l^ = - ^ [6.25 cos2x sin(2V' + 7r/3) - 8 sin2x cos(2V' + 7r/3) ] 
^GW = -I^ cos2x sin(2V> + 7r/3) 



These can be combined to give 
/o=^ (l + ^) sin(2^ + V3) 



/i = 3 ^1 + ^ 4 j ^^^^ cos(2'0 + 7r/3) — 6 cosx sin(2^ + 7r/3) 

/2 = 4\/3 sin 2x cos(2V' + it/3) - — f 25 + J cos 2x sin(2V' + 7r/3) 
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Figure captions 

Figure 1 

Normalised quadrupole moments J„ (in units of GMB?) for quadrupolar solar oscillations with surface 
displacement Cn(-R) = 1- 
Figure 2 

Signal function for oscillations in the Newtonian potential (a) m = 0, (b) m = ±1, (c) m = ±2. 
Figure 2 

Signal function for graviational waves produced by solar oscillations (a) m = 0, (b) m = ±1, (c) 

m = ±2. 

Figure 4 

Limits on detectable signal from LISA with signal to noise S/N=5. (a) Instrumental noise, (b) binary 
confusion noise (Bender and Hils 1997). 

Figure 5 

Solar background noise in velocity as determined from the GOLF experiment on SOHO. 
Figure 6 

Comparison of signal to noise (S/N) in LISA experiment for S/N=3 in solar velocity experiments, for 
oscillation modes with £ = 2, m = 2. 

Figure 7 

Minimum surface velocity amplitude for signal to be detectable by LISA with S/N=3 for modes with 
m = and m = 2. 
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